Test of exotic scalar and tensor interactions in K^,^ 
decay using stopped positive kaons 

A.S. Levchenko", A.N. Khotjantsev", A. P. Ivashkin", M. Abe^ 
M.A. Aliev'', V.V. Anisimovsky", M. Aoki^ I. Arai^ 
Y. Asano^, T. Baker'^''^, M. Blecher^, P. Depommier'^, 
^ : M. Hasinoff^ K. Horie'', H.C. Huang*, Y. Igarashi^ 

> ; T. Ikeda^ J. Imazato^ M.M. Khabibullin^ Y.G. Kudenko", 

2 : Y. Kuno^''^, L.S. Lee^ G.Y. Lim^'^ J.A. Macdonald^ 

D.R. Marlow^, C.R. Mindas"^, O.V. Mineev", C. Rangacharyulu^, 
S. Shimizu^, Y.-M. Shin^, A. Suzuki^ A. Watanabe^ 
N.V. Yershov", T. Yokoi^ 



O 



> 
00 

O 



Oh: 



(E246 KEK-PS COLLABORATION) 



"Institute for Nuclear Research of RAS, Moscow 117312, Russia; 
^University of Tsukuba, 305-0006, Japan; 
^KEK, 305-0801, Japan; 
'^University of Saskatchewan, S7N OWO, Canada; 
5J^! '^Virginia Polytechnic Institute and State University, USA; 

\ ■^University of Montreal, H3C 3J7, Canada; 

^University of British Columbia, Vancouver, V6T IZl Canada; 
^ ■ ''Osaka University, 560-0043, Japan; 

'National Taiwan University, Taipei, Taiwan; 
■'Institute of Physical and Chemical Research, Hirosawa 2-1, Japan; 
'^Korea University, Seoul 136-701, Korea; 
'TRIUMF, V6T 2A3, Canada; 
'Princeton University, NJ 08544, USA. 



rriT 



Abstract 

The form factors of the decay n^e'^u (-ft'es) have been determined from the 

comparison of the experimental and Monte Carlo Dalitz distributions containing 
about 10^ i^es events. The following values of the parameters were obtained: A+ = 
0.0278±0.0017(stat)±0.0015(syst), fs/f+{0) = 0.0040±0.0160(stat)±0.0067(si/st) 
and /t//+(0) = 0.019±0.080(stot) ±0.038(s?/st). Both scalar fs and tensor /t form 
factors are consistent with the Standard Model predictions of zero values. 
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1 Introduction 



The most general Lorentz invariant form of the matrix element of K^z decay can be 
written as Jl],!]: 

M (X /+(g')(Pi^ + P,o)^ne7A(l-75)^. 

+f^{q^)miUe{l - l^)uy + 2mx/5Me(l - 75)^1. 

+ {2fT/mK){PKf{P^oYueax^{l - 75)^., (1) 

where f±{q'^), fs and Jt are the vector, scalar and tensor form-factors, Pk and P-j^ 
are four-momenta of the and 7r°. The vector form-factors f±{q^) are assumed 
to be linearly dependent on the momentum transfer squared = {Pk — P-nY ^"^^ 
they can be represented by the equation f±{q^) = /±(0)(1 ± A±(g^/m^o))- Because 
of the small mass of the positron the part of the matrix element that depends on 
f-{q^) is negligible and there are just three free parameters of the theory: A+, fs 
and /t- Within the Standard Model (SM), due to VT-boson exchange, no terms 
other than those of a pure vector nature are expected. A possible contribution to 
fs and fx from electroweak radiative corrections is negligibly small and, therefore, 
nonzero values of fs and fx would signal new physics beyond the SM. Actually, 
such deviations from zero for these form factors were measured in []T],^, but were 
not confirmed in a recent experiment 0], where fs and fx were measured with 
sensitivity similar to Refs. |[T|,^. 

In this paper, we present results of a reanalysis of i^es data taken with the E246 
detector at the KEK 12-GeV proton synchrotron. The first result obtained using 
this data set was published in 

2 Experiment 

The experiment was performed using the set-up constructed to search for T-violation 
in 7r^fi~^u decay. The experimental arrangement is shown in Fig. |l|, and is 

described in detail elsewhere Kaons with Pk+ = 660 MeV/c are identified by 

a Cherenkov counter, slowed in a BeO degrader and then stopped in a target array of 
256 scintillating fibers located at the center of a 12-sector superconducting toroidal 
spectrometer. Charged particles from kaon decays in the target were tracked by 
means of multiwire proportional chambers at the entrance (C2) and exit (C3 and 
C4) of each magnet sector, along with the target and a scintillation ring hodoscope 
around the target. The momentum resolution of cXp = 2.6 MeV/c a.t p = 205 MeV/c 
was obtained using mono-energetic products from the two-body decay 7T~^7r^. 

The energies and angles of the photons from 7r° decays are measured by a CsI(Tl) 
photon detector consisting of 768 modules 0. The photon detector covers a solid 
angle of Sir steradians, with openings for the beam entry and exit and 12 holes for 
charged particles to pass into the magnet gaps. The photon energy was obtained by 
summing the cluster energy distributed among several crystals surrounding a central 
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Figure 1: The layout of the E246 set-up 



crystaL The position of the photon cluster is determined using an energy-weighted 
"center-of-gravity" method. To suppress accidental background from the beam, 
timing information from each crystal was used. A time resolution of 3.8 nsec (a) for 
the wide photon energy range of 10-220 MeV was achieved. The energy resolution 
was obtained to be aE /-E'7 — 3.0% at = 150 MeV. The invariant mass resolution 
of a^^ = 7.5 MeV/c^ at M^^ = 133.1 MeV/c^, and angular resolution of ao = 2.4° 
were obtained in the experiment @,||. 

The study of the i^es decay required a trigger condition different from that of 
the main experiment. The following trigger was used to accumulate K^^ events 

Ck X Fidi X T0F2i x 27 x GapVeto. (2) 

The Cherenkov condition Ck ensures that the beam particle is a K^. The coin- 
cidence between Ck, delayed by 2 nsec, and the fiducial counters (Fidi) for the 
charged particles eliminates the triggers due to decays in flight. The time of flight 
signal from the corresponding gap {T0F2i) completes the charged particle trigger 
condition. The trigger also required two hits in the Csl (27). The requirement 
GapVeto eliminates events where the charged particle may have lost energy in the 
Csl crystals around the i^'^ hole. 

The data were collected for two spectrometer field settings, B=0.65 and 
B=0.9 Tesla. 
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3 Analysis 



A'e3 events were selected using the following requirements. A clean hit pattern in 
the target and a delayed decay at least 2 nsec after the arrival time measured 
by Cherenkov counter suppressed the decays-in-flight. Cuts on the charged- 
particle momentum, p < 190 MeV/c, and the opening angle between and vr", 
6'e+7rO < 170°, effectively removed the Kt^2 events. Tr"*" decays-in-fiight were rejected 
by a cut on the track reconstruction < 14. The photon conversion events were 
suppressed by requiring single hits in the ring and fiducial counters and a single track 
in the target. Events with more than two hits in the electromagnetic calorimeter 
were rejected. The cut on the invariant mass was 75 < M^^ < 140 MeV/c^. The 
main criterion separating Ke3 positrons from K^s muons was derived from the time- 
of-flight (TOF) measurement between the fiducial counters and scintillating counters 
(T0F2) located at the exit of each spectrometer gap (see Fig. |l]). The mass squared 
of the charged particle was determined as 



where L and r are the length and time of flight of the charged particle between 
fiducial and T0F2 counters, and c is the speed of light in vacuum. The TOF 
resolution is (Jtof = 300 ps, and the mass squared spectrum is shown in Fig. ||. 
To separate positrons from muons, a cut on the mass squared —4500 < M^Qp < 
4500 MeVVc^ was chosen. 

Compared to the analysis of 0], improvements in reconstruction of the charged 
particle track were made. In the new analysis, information about the kaon stopping 
position from the target (x, y) and ring counter (z) was included in the momentum 
reconstruction routine. The reconstruction of the vr" kinematics was also improved 
by using tighter time windows in the Csl and improved calibration. The Monte 
Carlo routines properly included the correct kaon stopping distribution in the tar- 
get and the electromagnetic shower leakage effects in the Csl. This allowed us to use 
all Csl crystals for reconstruction of the 7r°, while in Ref. [Q events were rejected 
where either photon hit a crystal adjacent to the 12 charged-partilcle holes. Overall 
optimization resulted in approximately a factor of four increase in acceptance. For 
0.65-T magnetic field we extracted 102k good i^es events with background contam- 
ination of 0.21% from K^^s and 0.34% from Kt^2 decays. 

The distributions of the — 7r° opening angle and the opening angle between 
photons from vr^ '-ff for K^s decay is shown in Fig. ^. Energy spectra of the K^s 
decay are shown in Fig. ^. 

The extraction method for the parameters is based on comparison of the exper- 
imental and Monte Carlo Dalitz distributions. If one allows for the existence of the 
exotic interactions, the Dalitz plot density is given by [p!0[] : 



= P\{L/t-cY -1) 



(3) 



(4) 



where 
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Figure 2: (a) Momentum spectrum of from K^s decay; (b) /3 spectra of e'^ 
(/3 = :^); (c) reconstructed TOF mass squared of e"*"; (d) energy deposit of in 
T0F2 counter. Solid line corresponds to experimental and dashed line to GEANT- 
simulated events. 
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Figure 4: K^z decays: a) energy of 71 (-E^^ > -£^72)) b) energy of 72; c) invariant 
mass M^^; d) energy of 7r° (i?^o = E^^ + E^^). Solid line - experiment, dashed line - 
Monte Carlo. 
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B={ml){E^-E'j2) 

C = {mDE'jA 

El = (m^ + ml- ml)/{2mK) - E^ 

2^ _ {2mK/me)Rs + {me/rriK + 2{E^ - Ee)/me)RT 
(1 + {mJmK)RT) 

Rs = fsl f+ Rt = fr/ f+ 



The parameters of the decay can be extracted by selecting the minimum of the 
variable defined as 

X' = 2. J2iNr - ND + Nr ■ logi-f^), (5) 
i=i ' 

where n is the number of bins over the Dalitz plot, the number of Monte Carlo 
events in each bin, and N^^^ the number of experimental events in each bin. 

Within the Standard Model the only parameter to be obtained is A+. It can be 
extracted in a model-independent way from the dependence of /+. Using eq. |^, 
a/ p{Ee, -Ejr) oc /+ = /+(0)(1 + A+g^/m^o), and, therefore, A+ can be extracted from 
the ratio 

(iVe.p/iVmc(A+ = 0))i/2 = (1 + A+gVm^o). (6) 

The dependence of this ratio shown in Fig. ^ allows us to determine the value 
A+ = 0.0278 ±0.0016(stot). 

For extraction of all three parameters A+, fs and Jt, the Dalitz distribution 
fit was performed using the Pe+, 0.^0^+ variables. The angle between the pion and 
positron, ^7rOe+5 was preferred instead of pion energy E^^o in order to reduce a sys- 
tematic error related to the energy leakage in the electromagnetic calorimeter. The 
bin sizes chosen were APe+ = 3.125 MeV/c and A6'e+^o = 4.2°. Radiative correc- 



tions to the Dalitz plot were taken into account according to Ginsberg The 



between experimental and Monte Carlo Dalitz distributions was minimized by 



a program based on MINUIT [|T2[. The values obtained for A+, and the scalar and 
tensor form factors are 

A+ = 0.0278 ± 0.0017(stat) 
fs/f+{0) = 0.004 ± 0.016(stat) 

fr/f+iO) = 0.019 ± 0.080(stat). (7) 

Since our analysis relies on the proper detector responses used by the Monte Carlo 
simulation, comparisons and checks were made wherever possible. For example, to 
check the values obtained and to see possible unknown systematic effects, a similar 
analysis of the scatter plot 6-y-y vs. Pe+ was also performed. From this approach it 
was found that A+ = 0.0281 ± 0mi7{stat), /s//+(0) = 0.001 ± 0.018(stot), and 
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/t//+(0) = 0.007 ± O.lOO(stat). No significant variations of the A+ value were seen, 
and values of the form factors were consistent with zero within a 1 a uncertainty. 

The main sources of systematic errors were related to detector inefficiencies, mis- 
alignments of the detector elements, background contamination, and uncertainties 
connected to the Monte Carlo simulation. All these sources were studied and the 
estimations of the systematic errors are presented in Table |1|. 

4 Result 

The result presented here is based on about 10^ good K^.^ events. We have obtained 

A+ = 0.0278 ± 0.0017(stat) ± 0.0015(s?/st) 
fs = 0.0040 ± 0.0160(stat) ± 0m67{syst) 

fx = 0.019 ± 0.080{stat) ± 0.038{syst). (8) 

Using only one data set at B=0.65 T, we have improved statistical errors by a factor 
of 1.5 and systematic errors for A4., fs and /r were reduced by the factors 2.0, 2.1, 
and 2.4, respectively, compared to our previous result Q. This result is in agreement 
with the Standard Model prediction and there is no evidence for a deviation from 
zero for the values of scalar and tensor form factors. We expect to further improve 
our accuracy after completing the analysis of the Ke3 data accumulated at B=0.9 T. 



9 



Table 1: Systematic errors. 



Source A+ fs fj 



MWPC's misalignment 





0006 





0010 





0022 


MWPC's spatial resolution 





0002 





0010 





0010 


e"*" identification 





0007 





0022 





0004 


Bremsstrahlung 





0010 





0015 





0280 


cut point 





0001 





0010 





0018 


Csl barrel misalignment 





0005 





0042 





0240 


and My^ cut point 





0003 





0020 





0027 


7 — > e'^e" conversion 





0001 





0018 





0100 


Pile-up in the Csl 





0004 





0013 





0015 


Total 





0015 





0067 





0380 
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